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*  Specific  Contact  resistivity  p,  (Q-cm2),  which  can  be  extracted  from  measured  contact  resis¬ 
tance  Rf  (rf),  is  a  commonly  used  measure  of  metal-semiconductor  contact  quality.  It  is  indepen¬ 
dent  of  current  flow  and  contact  geometry,  and  depends  only  on  th^- transport  properties  of  the 
metal  to  semiconductor  junction.  Unfortunately,  extraction  of  ^obtained  from  test  structures 
differing  in  size  or  design  often  disagree  by  as  much  as  an  order  of  magnitude  primarily  because 
simple  1-D  models  which  are  used  to  extract  the  value  of  pe  can  not  account  for  fringing  resistance 
associated  with  the  2-D  nature  of  the  current  distribution  around  the  contact  Jlf.  ‘  For  example, 
in  a  widely  used  model  [2]  for  cross  bridge  Kelvin  resistor  the  expression  pe  =  (Re  x  Contact  area) 
is  used.  But  in  practice,  a  sublinear  behavior  is  often  observed  in  a  plot  of  pc  vs  contact  area,  r 
from  which  a  meaningful  value  of  pc  can  not  be  extracted  if  this  1-D  model  is  used,  (3). 

In  this  work  we  have  used  a  2-D  simulator  to  analyze  three  commonly  used  contact  resistivity 
measurement  devices:  the  cross  bridge  Kelvin  resistor,  the  transmission  line  tap  resistor,  and  the 
contact  end  resistor,  which  are  all  shown  in  Fig.  1.  This  analysis  was  used  to  assess  the  extent  of 
the  error  introduced  by  the  use  of  1-D  models  for  extraction  of  pe,  and  to  develop  a  method  by 
which  a  2-D  model  may  be  used  to  extract  a  more  reliable  value.  This  method  is  then  generalized 
so  that  a  limited  number  of  simulations  may  be  used  to  extract  pe  from  devices  with  a  wide  variety 
of  geometries  and  sheet  resistances.  The  technique  has  been  used  to  obtain  pe  for  PtSi,  Pd2Si,  W 
and  A1  contacts  to  N+  and  P+  Si  with  dopant  concentration  in  the  range  of  101S  -  1020  cm-3. 

Fig.  2  shows  the  results  of  these  simulations,  where  the  Helmholz  equation  y2V=V/(p(./ R,) 
was  solved  in  the  contact  region,  and  the  Laplace  equation  was  solved  in  the  remainder  of  the 
diffusion  region.  For  the  ideal  structures  where  the  contact  width  is  the  same  as  that  of  diffusion 
(6=0)  the  conventional  1-D  models  shown  in  Fig.  1  give  the  same  results  as  the  2-D  models. 
However,  for  real  structures  with  finite  values  of  6  the  measured  values  of  the  contact  resistance 
(Rjfc,  R/  and  Re  for  Kelvin,  transmission  line  and  end  resistances,  respectively)  are  significantly 
larger  than  the  1-D  model  prediction.  This  is  due  to  the  added  voltage  which  is  generated  by 
current  flowing  in  in  the  overlap  region  between  the  contact  and  the  diffusion  edge.  This  effect 
becomes  larger  for  lower  pe,  higher  6  and  higher  sheet  resistance.  Clearly  the  use  of  1-D  models 
will  result  in  overestimation  of  pf.  This  effect  is  most  serious  for  the  contact  end  resistance 
structure,  and  least  serious  for  the  transmission  line  tap  structure. 

By  systematically  varying  either  the  diffusion  width  or  the  contact  area  the  fringing  parasitic 
resistance  can  be  isolated  from  the  true  interface  contact  resistance  by  comparing  the  experimental 
data  to  2-D  simulations  as  shown  in  Fig.  3  for  the  cases  of  Kelvin  and  contact  end  resistance 
structures.  By  simulating  the  same  device  over  a  variety  of  resistivities,  and  by  fitting  the  data 
points  to  one  of  the  simulated  curves,  pt  can  be  determined.  We  have  fabricated  the  three 
structures  on  the  same  test  chip  and  demonstrated  that  the  same  value  of  pe  is  obtained. 

Unfortunately,  this  method  requires  a  new  set  of  simulations  for  each  value  of  pe  and  R,.  But 
by  using  the  previously  developed  theory  of  scaling  of  contacts  [1,4]  we  have  determined  that  a 
more  efficient  extraction  technique  is  to  express  the  ratio  Rr/R,  (where  Rf  is  actually  Rfc,  R/, 
or  R,)  as  a  single  dimensionless  parameter.  This  parameter  is  plotted  against  any  appropriate 
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geometrical  variable,  such  as  contact  length  or  contact  area,  resulting  in  a  set  of  universal  curves 
which  apply  to  a  very  wide  range  of  pe  and  R,  values.  This  method  is  illustrated  in  Fig.  4.  and 
is  applicable  even  to  the  “chain  of  contacts”  structure,  since  that  is  simply  a  special  case  of  the 
transmission  line.  If  contact  length  and  the  parameter  Lt  are  normalized  to  the  overlap  width  <5, 
then  these  general  curves  will  also  apply  to  a  wide  range  of  contact  and  diffusion  sizes. 

The  curves  shown  in  Fig.  2  reveal  that  of  the  three  test  devices  simulated,  the  end  contact 
resistor  is  the  most  sensitive  to  diffusion  overlap  6.  Varying  this  overlap  from  0  to  2.5  microns  can 
cause  the  contact  resistance  to  vary  by  two  orders  of  magnitude  or  more.  Therefore  this  device 
cannot  be  relied  upon  to  provide  accurate  estimates  of  pe.  Of  the  two  remaining  structures,  the 
transmission  line  is  the  least  sensitive  to  variations  in  geometry,  but  in  practice  it  can  be  inac¬ 
curate,  because  of  the  diffusion  leading  up  to  the  contact.  Because  this  diffusion  is  part  of  the 
current  path,  and  also  part  of  the  voltage  tap,  its  resistance  must  be  subtracted  by  extrapolation 
from  the  total  measured  resistance.  This  diffusion  resistance  can  typically  be  2  orders  of  magni¬ 
tude  larger  than  R/,  rendering  the  subtraction  very  unreliable.  We  would  therefore  recommend 
the  cross-bridge  Kelvin  resistor  as  the  most  accurate  of  these  structures,  because  it  has  a  mini¬ 
mum  parasitic  resistance  in  its  voltage  tap,  and  is  therefore  relatively  insensitive  to  variations  in 
geometry.  When  used  in  conjunction  with  a  2-D  model,  it  can  provide  very  accurate  results. 
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Fig  1.  Test  structures  for  contact  resistance  measurements  along  with  1-D  model 
expressions:  (a)  transmission  line,  (b)  end  contact  and  (c)  cross  bridge  Kelvin 
resistance  structures. 
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Fig  2.  2-D  simulations  of  the  test  structures  showing  the 
impact  of  the  diffusion  collar  width,  6,  around  the  contact 
on  the  measured  contact  resistance,  for  (a)  Kelvin,  (b) 
end  contact  and  (c)  transmission  line  resistance  structures. 
The  1-D  model  is  identical  to  the  case  where  6  =  0.  The 
error  is  stronger  for  lower  pc,  higher  6  and  larger  contact 
size. 
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Fig  3.  Illustration  of  the  method  used  to  extract  pc  by  plotting  Rt  data  against 
geometry  and  fitting  the  data  to  the  nearest  simulation  curve  for  (a)  Kelvin  resis¬ 
tance  and  (b)  end  contact  resistance  structures.  The  dashed  line  in  part  a  shows 
the  resistance  given  by  the  1-D  model  for  Kelvin  resistance. 


Fig  4.  Generalized  technique  of  pc  extraction  for  the  Kelvin  resis¬ 
tor  structure  by  comparing  the  experimental  data  with  a  univer¬ 
sal  plot  olRjR.  vs  contact  area  as  a  function  of  transfer  length 
(£1  =  \/ Pc  l  R>)- 
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2-D  Model  of  Metal/Diffusion  Contact 


variation  in  depth  ignored 


system  described  by  potential  in  diffusion  ‘sheet',  V 
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Universal  Curves  from  Seal 
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